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I review the current status of lattice QCD calculations and the progress that we have made in 
consolidating the 'lattice QCD revolution' of four years ago. Significant results from formalisms other 
than the improved staggered formalism, which has been leading the revolution, are now appearing. 
Comparison between formalisms gives additional confidence in the results. 

More precision tests against experiment have been made and predictions that are needed by the 
experimental programme have been improved. There has been particular progress this year in charm 
and strange physics. 
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INTRODUCTION 

QCD is a key component of the Standard Model of par- 
ticle physics. On the one hand, it gives us a rich spectrum 
of bound states of quarks and giuons whose properties 
are predictable from QCD if we can solve the theory. On 
the other hand, the confinement of quarks comphcates 
the determination of the properties of quarks from ex- 
periment because only hadrons can be studied directly. 
Lattice QCD enables QCD effects to be calculated 'from 
first principles' in the hadronic regime where the theory 
is strongly-coupled and nonlinear. Accurate results can 
provide stringent tests of QCD when compared to exper- 
iment as well as providing input to our understanding of 
the Standard Model in the quark sector. 

In current lattice QCD we can calculate the masses of 
'gold-plated' hadrons (those with small width well below 
Zweig-allowed decay thresholds) and simple decay matrix 
elements that include at most one gold-plated hadron in 
the final state. There arc both statistical and system- 
atic errors from lattice calculations, however, and it is 
important to understand the sources of these so that an 
optimal strategy can be developed to minimise the total 
error from the calculation. 

One area in which lattice QCD is making an important 
contribution is that of flavor physics and the determina- 
tion of elements of the Cabibbo-Kobayashi-Maskawa ma- 
trix. This is linked to the worldwide experimental pro- 
gramme determining weak decay and mixing rates for 
bottom and strange hadrons with a view to fixing all 3 
sides of the unitarity triangle and so testing the internal 
consistency of the Standard Model's description of CP 
violation. The experimental programme will achieve er- 
rors of a few percent on the decay rates (and has already 
done this on mixing rates) and needs theoretical input for 
the Standard Model prediction to extract the appropriate 
CKM element. The final error for the CKM element, and 
corresponding side of the unitarity triangle, is currently 
limited by the error from lattice QCD and this must be 
reduced to be much closer to the experimental error, if 
that accuracy is not to be wasted. 

The good news is that it now looks possible to achieve 



this as a result of key advances over the past decade in 
understanding how to discretise QCD accurately onto a 
space-time lattice. Lattice QCD calculations are numer- 
ically extremely expensive, but the savings that have re- 
sulted from improved discretisation have at last brought 
realistic calculations within the power of current day su- 
percomputers [l|. I will outline the advances that have 
made the recent calculations possible and their implica- 
tions for future work, concentrating on calculating rel- 
evant to the flavor physics and CKM programme. The 
Proceedings of this year's lattice conference should be 
consulted for a more general view [2| . 



LATTICE CALCULATIONS 

Lattice QCD proceeds by the numerical evaluation of 
the Fcynman path integral [3[. For this integral to be 
finite and well-behaved we work with a finite volume 
of 4-dimcnsional space-time in which time is rotated to 
(imaginary) Euclidean time. The space-time within our 
volume is split up into a lattice of points with lattice 
spacing a and there are then a finite number of quark 
and gluon quantum fields residing on the sites (quarks) 
and links (giuons) of the lattice. We then have to in- 
tegrate over all possible values of these fields, weighted 
exponentially by (minus) the action, S (integral of the 
Lagrangian), of QCD. In practise this means using a ran- 
dom process to generate sets of possible gluon fields, one 
for every link of the lattice, called configurations. This is 
the 'data generation' phase of a lattice QCD calculation. 
If we generate these configurations with probability e~^ 
then we are preferentially choosing configurations that 
contribute most to the path integral and we can evaluate 
it efficiently. We call these configurations 'typical snap- 
shots of the QCD vacuum'. A set of configurations for a 
particular set of QCD parameters is called an ensemble. 
An ensemble for a good calculation will generally have 
several hundred configurations in it. The evaluation, or 
analysis, stage of a lattice QCD calculation consists of 
'measuring' various functions of the gluon fields that cor- 
respond to a particular observable, such as a correlation 
function from which a hadron mass can be determined. 



The function of the gluon fields is evaluated on each con- 
figuration of the ensemble and the mean value and its 
statistical error determined. The statistical error will 
depend on the number of configurations in the ensem- 
ble, i.e. how well the path integral is approximated by 
this procedure. The statistical error for a 'measurement' 
varies as the inverse square root of the number of config- 
urations in the ensemble and to reduce this error to 1% 
is perfectly feasible in current calculations for quantities 
such as flavour non-singlet ground-state hadron masses. 
Reducing systematic errors to this level is much harder, 
and it has taken many years of effort to understand how 
to do this. 

The key source of systematic error is that coming from 
the discretisation of the Lagrangian of QCD onto a space- 
time lattice. Discretisation errors invariably arise when 
equations from continuous space-time are discretised for 
numerical solution. Typically physical results will de- 
pend on the unphysical step-size or spacing chosen for 
the discretisation. The physical result for continuous 
space-time will be obtained either by extrapolating in 
the step-size to zero, or reducing these errors to a known 
and acceptable level. The errors arise, for example, from 
the approximation of derivatives by finite differences, and 
for classical equations they are readily reduced by using 
higher-order differencing schemes which correct for the 
errors being made in low-order schemes by adding addi- 
tional terms to cancel them. This has the effect of raising 
the power of the step-size at which errors first appear. 
This therefore reduces them, at fixed small step-size, or 
allows you to achieve you the same errors with a cheaper 
calculation with larger step-size. 

Exactly the same considerations apply in lattice QCD. 
The discretisation step-size is the separation between 
points in the lattice and is known as the lattice spac- 
ing, a. The number of lattice points required in a fixed 
physical volume grows as a~^ as the lattice spacing is 
reduced. However, the cost of a calculation grows much 
more rapidly than this, approximately like a~^, because 
the efficiency with which statistically independent con- 
figurations can be generated also falls as a function of a. 
This makes it impossible to work at very small values of 
a and hence understanding the systematic discretisation 
errors associated with the lattice spacing has been criti- 
cal in enabling us to obtain results at values of a at which 
we can afford to do the calculation. 

Discretisation of the gluon piece of the QCD action, in- 
cluding the effect of 'improvement' to reduce its discreti- 
sation errors is relatively straightforward. It is the quark 
piece of the QCD action that causes the most headaches 
and controversy. Since quark fields anticommute, as they 
are fermions, they cannot be included directly in our 
numerical simulations with commuting arithmetic, but 
must be integrated out of the Feynman path integral. 
The path integral then becomes an integral over gluon 
fields only, with the effect of quarks being implicit as 



functions of the gluon field. If we write the quark piece 
of the QCD Lagrangian as 



Lq,QCD = i'{l ■ D + m)^ = iiMij} 



(1) 



then we call M the 'fermion matrix'. 7 are the Dirac 7 
matrices, D is the covariant derivative that includes cou- 
pling to the gluon field and m is the quark mass. The 
quark field, "0, is a 4-spin, 3-color vector on every site 
of the lattice so M is a 121/ x \2V matrix where V is 
the number of sites on the lattice. The result of inte- 
grating quarks out of the Feynman path integral means 
that their effects appear in two different ways. For va- 
lence quarks we must calculate the 'quark propagator', 
M~^, and combine quark propagators to make hadron 
correlation functions. For sea quarks, we need to include 
dct(M) in the generation of gluon field configurations to 
make 'typical snapshots of the vacuum' that include the 
effect of sea quark-antiquark pairs being produced by en- 
ergy fluctuations in the vacuum. 

The calculation of rows of M~^ is numerically costly 
because M is such a large matrix. It is also ill-conditioned 
as TO ^ and this is unfortunate since we have quarks 
in nature, the u and d quarks, that are very light. The 
inclusion of det(M) in the generation of gluon field config- 
urations is numerically even more costly since it involves 
many calculations of rows of M~^. Here the fact that the 
u and d quarks are so light (and s to a lesser extent) is 
even more of a problem, but at the same time it is exactly 
the reason that they are the most physically important 
as sea quarks since they can readily be generated by a 
vacuum energy fluctuation. 

The quarks that we need to be concerned with are only 
the 5 lightest, because the t quark does not form hadronic 
bound states, c and h quarks are themselves too heavy 
to have any significant effect as sea quarks, and appear 
in existing calculations only as valence quarks, u, d, and 
s sea quarks are physically important, however, as stated 
above. Their inclusion is necessary, for example, for the 
strong coupling constant to run correctly. Without this, 
quantities which are sensitive to different energy scales 
will not agree. Early lattice calculations did not have 
the computer power to include sea quarks and so they 
were dropped in what was known as the 'Quenched Ap- 
proximation'. This approximation has systematic errors 
at the 10% level which can be moved around between 
different quantities depending on how the calculation is 
done, because it is not an internally consistent approxi- 
mation to QCD. Subsequently calculations including sea 
quarks were done, but they included only two flavours 
of sea quarks (i.e. u and d) but with masses that were 
10-20 times too big (i.e. around the mass of the strange 
quark). The numerical cost of including light u and d 
quarks means that it is very important to find a fast and 
accurate discretisation of the quark Lagrangian so that 
affordable calculations can be done at moderate values 
of a (around O.lfm). In practice calculations are still 



done at multiple values of the u and d quark masses that 
are too heavy and then extrapolations are made to the 
physical point. These extrapolations can be guided by 
chiral perturbation theory [3|, which gives an expansion 
in powers of the u/d quark masses, provided that we are 
close enough to 77i„/d = 0. 

Controversy enters at this point since there are many 
different formulations for quarks in lattice QCD, depend- 
ing on the approach taken to solving the infamous 'dou- 
bling problem'. When the Dirac Lagrangian above is dis- 
cretised onto the lattice it describes 2'^ continuum quarks 
in d dimensions rather than 1. This gives 16 quarks in 4 
dimensions. We call the 15 additional quarks 'doublcrs' 
or additional 'tastes' of quark. Different formulations 
take a different attitude to this problem with consequent 
effects on speed and discretisation errors. Maintaining 
as much as possible of the chiral symmetry of QCD, un- 
der which left- and right-handed projections of the quark 
field can be separately rotated in the absence of a quark 
mass, is also important. The pion mass is guaranteed 
to vanish at zero quark mass, for example, because it 
is the Goldstone boson of spontaneously broken chiral 
symmetry. If this can be maintained on the lattice it 
is a big advantage in terms of easily being able to find 
your way, from a given lattice quark mass, to the phys- 
ical point for u and d quarks, close to zero quark mass. 
In some lattice quark formalisms this property has to be 
given up and then the value of the quark mass (which 
can be negative) corresponding to zero pion mass has 
to be searched for, which is an added complication. On 
the other hand maintaining the best possible version of 
continuum chiral symmetry on the lattice is numerically 
extremely expensive and not a particular advantage for 
a lot of calculations. 

Tabic [J gives a brief outline of the existing formalisms 
for which significant numbers of gluon field configura- 
tions including the effect of sea quarks have been made, 
and a significant amount of physics analysis done. It in- 
cludes the names of the collaborations chiefly using those 
formalisms. It is also possible to use a different valence 
quark formalism from that used for the sea quarks in- 
cluded in the gluon field configurations that you are work- 
ing on. For example, several collaborations have used 
domain wall valence quarks on the MILC collaboration 
configurations with improved staggered sea quarks. The 
table briefly describes the good and bad points of each 
formalism with respect to speed and chiral symmetry. 
There are other technical issues associated with each for- 
malism that the collaborations involved have spent con- 
siderable time understanding. In most cases they amount 
to restrictions on how the continuum (a —^ 0) and chi- 
ral (mq — !■ 0) limits are approached, and they limit how 
small a value can be taken for m^ at finite a. For exam- 
ple, for improved staggered quarks no attempt is made 
to solve the doubling problem and the doublcrs are taken 
care of by effectively 'dividing by 16'. Away from the 



a = limit that this amounts to a discretisation error 
that can cause problems if the valence quark mass is 
taken too light. This has been demonstrated both numer- 
ically and in 'staggered chiral perturbation theory', which 
takes account of the effect of the doublcrs in chiral pertur- 
bation theory. For a recent review of issues for staggered 
quarks see [J] . The discretisation errors are reduced fur- 
ther in a new highly improved staggered quark (HISQ) 
formalism introduced recently by the HPQCD collabora- 
tion [51] , that will be discussed further below. For domain 
wall quarks a small additive quark mass renormalisation 
appears which means that care must again be taken with 
small valence quark masses. For a recent review of is- 
sues for domain wall quarks see [61] and for the related 
overlap quarks see [7|. For twisted mass, see [8[ and for 
clover [9|. In practice these issues are under control in 
existing calculations. Another important problem, and 
one on which there has been a lot of recent work, is that 
of the lattice volume. 

The dependence on rriu/d hi chiral perturbation is ob- 
tained by integrating over the momenta of virtual pious. 
On the lattice the momenta are restricted to the dis- 
crete set available on a particular spatial volume. The 
lowest momentum (from discretising a derivative as a 
simple symmetric finite difference) in units of the lat- 
tice spacing is given by pa = sin(27r/L) where L is the 
number of lattice points on a side of the lattice. Thus 
the smallest momentum (and therefore the infrared cut- 
off on the integral) is inversely proportional to the lattice 
size. At large quark masses the integral is cut-off by the 
quark mass but at small quark masses it becomes sen- 
sitive to the volume and this will distort the behaviour 
as a function of quark mass. This means that care must 
be taken at small u/d quark masses to be working on a 
large enough spatial volume (that can be estimated using 
finite volume chiral perturbation theory ll| ) , or to have 
multiple volumes and check the volume dependence ex- 
plicitly (sec, for recent examples [sl, |lO[). Larger spatial 
volumes are more expensive numerically, of course, since 
at fixed lattice spacing they require more lattice points. 

Figure [T] shows the status of configurations that exist 
and are being generated as reported at the Lattice 2007 
meeting [2|. Improved staggered quarks have the best 
coverage of different values of the lattice spacing, but 
other formalisms are making good progress in generating 
configurations with light u/d masses on relatively fine 
configurations. Not all of these include sea s quarks at 
present but all have plans to do so. 

From the very brief discussion above it will be clear 
that several different formalisms to handle quarks on 
the lattice are possible and different collaborations have 
strong views on the optimal approach. Of course, all the 
formalisms should give the same physical result in the 
end, and agreement between different formalisms is valu- 
able confirmation of the validity of lattice QCD. Different 
formalisms have particular strengths and so it tends to be 
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TABLE I: Different quark formalisms currently being used 
for generating gluon field configurations including the effect of 
sea quarks. Collaborations currently using these formalisms 
are given on the right. 
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FIG. 1: Configurations made by various collaborations using 
different quark formalisms and including the effects of u/d sea 
quarks only (n/ = 2) or u,d and s (n/ = 2-|- 1). The axes are 
chosen so that is clear how close to the 'real world' (circled at 
a = and m-„ = 0) the configurations are. The x axis gives 
a? since systematic errors for the actions shown here appear 
first at that order in a. The y axis gives m% for n mesons 
made of quarks with the same quark mass as those in the sea 
for the ensemble with lightest sea quark mass at that value 
of a. This is a measure of how far from the chiral limit the 
configurations are since m^ is proportional to the light quark 
mass. Other considerations for the quality of an ensemble are 
the spatial volume of the lattices and how many independent 
configurations there are. Open symbols give parameters for 
configurations currently being generated. 



true that collaborations using different formalisms have 
a different focus for their physics programmes. 

There is overlap, however, in a lot of very basic calcula- 
tions that act as a test case for the results. One example 
is the calculation of the pion decay constant /^ that pa- 
rameterizes the purely leptonic decay of a charged tt to 
leptons via a W boson (sec figured]). The leptonic decay 



rate is proportional to the square of /tt multiplied by the 
square of Vud which is the appropriate CKM element by 
which the u and d in the 7r+ couple to the W . Given 
a very accurate value for V^d from superallowed (3 de- 
cay it is possible to determine /tt from the experimental 
leptonic decay rate and compare lattice QCD results to 
this. The calculation of /tt is relatively simple in lattice 
QCD for formalisms that have enough chiral symmetry 
to have a partially conserved axial current (so that no 
renormalisation is required). Then the calculation of /^ 
is straightforwardly obtained from the same correlators 
that are used to obtain niT, . Calculations need to be done 
at several different values of the u/d quark mass and the 
lattice spacing a. For light enough rriu/d ^ fit of the 
results to chiral perturbation theory can be done that al- 
lows an extrapolation to the physical value of m.^^u (that 
gives the physical value of tti.^. ) . The different values of a 
allow for a simultaneous extrapolation to the continuum 
limit. At the physical point a comparison can be made 
with the value obtained from the experimental rate. Us- 
ing improved staggered quarks or HISQ agreement with 
experiment is obtained with 1.5% errors [12| . This year 
new results from the twisted mass collaboration also gave 
1% accurate /tt values (without fixing the lattice spacing) 
in the continuum and chiral limits with 2 fiavors of sea 
quarks at a range of m„/(i and with two values of the 
lattice spacing [13| . 




FIG. 2: Annihilation of a charged pseudoscalar meson to a 
W boson. The rate for this can be measured experimentally 
and depends on the appropriate CKM element coupling the 
annihilating valence quark and antiquark and on the 'proba- 
bility for the quark and antiquark to be in a position inside 
the meson to annihilate'. This probability is square of the 
matrix element of the axial vector current between the meson 
and vacuum, parameterised by /.^. It must be calculated in 
lattice QCD because it is sensitive to the long distance QCD 
effects that confine the quarks in the meson. 

A comparison of different quark formalisms can also 
be made away from the chiral and continuum limits if 
discretisation errors are small. This comparison for for- 
malisms with good chiral symmetry is shown in Figure [3] 
and shows very encouraging agreement. To make a more 
detailed comparison would require more work to remove 
effects of finite volume that different groups have han- 
dled in different ways. /,r is sensitive to finite volume 
effects at the lower end of the quark masses shown here 
and this has been explored by explicit calculations for 



the MILC [101, ETMC H and RBC/UKQCD [ij results 
shown here, and compared to chiral perturbation theory 
expectations. The conclusion is that the effects are of the 
size you would expect - for example a 2% shift downwards 
of /,r on a 2.5fm lattice size (compared to the infinite vol- 
ume real world) for m^ « 250 MeV. In FigureOthe MILC 
lattices at that 77i,r are 2.8fm, the RBC/UKQCD lattices 
are 2.7fm (using their 24'^ lattices rather than earlier 16'^ 
ones) and the ETMC lattices are 2.2fm. The NPLQCD 
results use the MILC coarse lattices with domain wall 
valence quarks |15| . Results with the clover formalism 



are also available but do not at present seem to ag ree as 
well, perhaps because of renormalisation issues [16[ . It is 
very exciting that we are now able to do a comparison 
between formalisms at this level of precision and more of 
this will become possible in future years. 
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FIG. 3; A comparison of results for the pion decay con- 
stant obtained from different quark formalisms in lattice QCD 
including the effect of sea quarks. The results are plotted 
against m^ (which is proportional to m^/ci) and each calcula- 
tion needs results at multiple values of a and m^i^ so that an 
extrapolation to the physical point where m^/^ is small and 
a = can be made. Some formalisms have results only at one 
value of a so far. The lines give the chiral extrapolation for 
the results from using the highly improved staggered quark 
formalism (HISQ) on the MILC ensembles using improved 
staggered quarks at each of 3 values of a plus the a = line 
obtained [l2|- The results show very encouraging agreement 
between the different formalisms. The experimental point 
marked uses the leptonic decay rate of the n meson and Kid 
from the particle data tables obtained from super-allowed /3 
decay. 

Handling heavy quarks on the lattice raises rather dif- 
ferent issues from that for light quarks. If we use one 
of the standard light quark formalisms discretisation er- 
rors will be set by powers of viqa where toq is the heavy 
quark mass, rather than, as is typical for light quark 
quantities, powers of Aqcdo.- For the b quark mta » 1 
on typical lattices and so no amount of improvement 



can control the discretisation errors. However, we can 
make use of the fact that the b quark is nonrelativistic 
inside its bound states and that nib is just an overall 
mass scale that does not affect the internal dynamics 
very much. The HPQCD collaboration has done a lot 
of work on bottomonium and B physics using the non- 
relativistic effective theory called NRQCD discretised on 
the lattice |17| . An alternative is to expand in powers 



of l/rriQ away from the infinite quark mass (static) limit 
in the HQET approach to B physics. The FNAL/MILC 
collaboration use a rclativistic formalism (the clover for- 
malism) making use of nonrelativistic understanding to 
remove key discretisation effects [l8|. Again all the dif- 
ferent formalisms should give the same physical results. 
For c quarks the situation is less clear because on typical 
lattices rrica « 0.5. The FNAL/MILC collaboration have 
made strong use of their clover 'Fermilab' formalism for 
charmonium and D physics. The HISQ formalism uses 
a variant of the improved staggered formalism to reduce 
further the discretisation errors associated with multiple 
tastes. Other more standard discretisation are already 
removed to a very high level for this action and so it gives 
excellent results for c physics that will be discussed be- 
low. It can also, of course, be used for m, d and s quarks, 
both as valence quarks (see below for results) and in the 
sea (in progress) as a further check that lattice discreti- 
sation errors are well understood. 



LATTICE RESULTS 2007 

As discussed in the introduction lattice QCD has an 
important contribution to make to determining the el- 
ements of the CKM matrix. For each CKM element 
there is a gold-plated clcctroweak decay or mixing pro- 
cess whose rate, as for it leptonic decay, will be (up to 
known kinematic factors) the product of that V^f^ and 
the square of a lattice QCD amplitude given as a decay 
constant, form factor or bag parameter that expresses 
the probability of the quarks confined inside the meson 
undergoing that process. The CKM matrix is given in 
Figure 2] with the corresponding leptonic, semileptonic 
and mixing processes for each element. 

Of course, it is not sufficient to calculate only these pro- 
cesses in lattice QCD. It is important to have a number 
of cross-checks against other processes that are similar 
and well-known experimentally, for example electromag- 
netic decay rates, as well as checking a variety of hadron 
masses. So a complete programme of this kind encom- 
passes the whole range of flavor physics. Figure [5] shows 
a 2007 update of a range of quantities obtained from 
lattice QCD calculations with improved staggered sea 
quarks [l|. The impressive agreement across the board 
provides strong confirmation that lattice QCD is accu- 
rately describing the real world when sea quarks arc in- 
cluded. A companion plot of results in the quenched ap- 
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FIG. 4: The CKM matrix with corresponding gold-plated 
processes that allow the value of each element to be deter- 
mined by combining experiment with a lattice QCD calcula- 
tion. 



proximation showing 10% errors has now been dropped 
since there is no longer any point in attempting to pro- 
duce quenched results for comparison. It is to be hoped 
that other formalisms will produce their version of this 
'ratio plot' soon. 
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FIG. 5: The points show lattice QCD results divided by 
experiment for a range of quantities from light quark physics 
to bottomonium physics, compared to the right answer of 1.0. 
f-n and JK are the tt and K decay constants, described further 
in the text; mjv is the nucleon mass. T(1P — IS) denotes, 
for example, the difference in mass between the lowest Xb 
states (spin- averaged) and the T. ^{IP — IS) is the same 
quantity for charmonium. The results come from analysis by 
the FNAL, HPQCD and MILC collaborations on the MILC 
ensembles that include improved staggered sea u, d and s 
quarks 
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The way in which lattice QCD calculations are done 
was briefly described earlier. Here it is useful to describe 
further how the parameters of QCD are fixed. QCD de- 
scribes an enormous range of physics with very few pa- 
rameters: a mass for each quark flavor and a coupling 
constant. For current lattice QCD calculations the t 
quark is ignored and the u and d quarks arc taken to have 



the same mass for numerical speed and convenience. This 
gives 5 parameters to be fixed and we do this by fixing 5 
hadron masses or mass differences to their experimental 
values. It is important to use gold-plated hadrons since 
hadrons that decay strongly or are close to decay thresh- 
olds will be sensitive to coupling to real or virtual decay 
channels that will distort the mass and, if there are sys- 
tematic errors here, these will then be fed into the rest 
of the calculation. The hadron mass being used should 
be sensitive to the quark mass it is being used to fix but 
preferably not sensitive to other quark masses to avoid a 
complicated iterative tuning problem. For the improved 
staggered results shown here the lattice spacing is fixed 
from the radial excitation energy in the T system, i.e. 
the difference in mass between the T' and the T, which 
turns out to be insensitive to all quark masses [19| . The 
u/ d quark mass is then fixed from Af^, the s quark mass 
from Mk [IJ, l20|, the c quark mass (using HISQ) from 
Mrj^ \\% and the b quark mass (using NRQCD) from 
Mx [19[ . Other gold-plated quantities can then be calcu- 
lated with no free parameters and these are shown in Fig- 
ure[51 Other choices to fix the parameters could be made, 
particularly for the lattice spacing itself. The ETMC col- 
laboration advocate using /^ [8| and the RBC/UKQCD 
collaboration mo [6[. At intermediate points in the cal- 
culation it is convenient to use a nonphysical quantity 
to determine the relative lattice spacing very accurately 
between different ensembles. Most groups use a distance 
parameter from the heavy quark potential, either tq or 
Ti corresponding to different values for the force between 
two infinitely massive quarks [23[. This can be deter- 
mined with better than 0.5% statistical accuracy but its 
physical value cannot be directly determined from exper- 
iment so at the end there must be conversion to physical 
units using, for example, T{2S — \S) [l9|. 

Lattice QCD then provides a very natural and accurate 
way to determine the parameters of QCD, superior to any 
other method, and results from this have made their way 
into the particle data tables [27| . Further work on this is 
ongoing, but I will not report on it here. There will be a 
number of new results next year with improved accuracy 
for quark masses. 

Here I will concentrate on lattice QCD calculations 
for CKM element determination starting with Vus- The 
determination of f-„ was described above and there is 
an analogous calculation for the K meson, yielding fx- 
Again an 'experimental' result for Jk can be obtained 
from the experimental leptonic decay rate combined with 
a Vus from elsewhere. Usually Vus is taken from K 
semileptonic decay which I will discuss shortly. This is 
the experimental result that is used in the ratio plot of 
Figure [3 (Note that the experimental value for Jk used 
there has been updated from that quoted in the 2006 par- 
ticle data tables [27| to be consistent with their quoted 
value of Kis ) • The lattice results in that plot come from 
HISQ valence u/d and s quarks on the MILC improved 



staggered ensembles, with lattice errors of 1-2%. The 
calculation of Jk/ fir in lattice QCD can be done with 
a smaller error - 0.6% for HISQ on the MILC ensem- 
bles - and this can be used, along with the ratio of the 
experimental leptonic deca y r ates 24 l27l | , to determine 



Vus/Vud and therefore Vus |25[. In this way the HPQCD 
collaboration recently obtained Vus = 0.2262(14) [I3 and 
the MILC collaboration updated their previous fx/ f-i^ 
analysis ^ to give 0.2246(-h25-13) ^. Both are com- 
petitive with the result quoted from Kli decay in the 
particle data tables of 0.2257(21). Figure [S] shows a com- 
parison of Vus values obtained from recent lattice results 
for JkI fir using 2-1-1 flavors of sea quarks. It is clear that 
this is an accurate way of determining Vus ■ The final re- 
sult is still completely dominated by theoretical error, 
however. To improve it further will require a number of 
improvements to the lattice calculation, chiefly working 
on larger volumes to reduce the finite volume error in /^r 
and reducing the uncertainty in the determination of the 
lattice spacing. 



RBC/UKQCD dw 
PACS-CS clover 

HPQCD HISQ 

on asqtad sea 

MILC asqtad 



0.22 0.225 0.23 

FIG. 6: Determinations of Vus from lattice QCD results for 
Jk / f-n combined with experimental results for the ratio of 
leptonic decay rates. The lattice QCD results use different 
formalisms [3, [I^, [H, H but all with 24-1 flavors (i.e. u/d 
and s) of sea quarks. 



Vus can also be determined from K semileptonic decay 
to ttIv. Because this is now a 3-body decay the quantity 
that must be determined in the lattice QCD calculation 
is a form factor that depends on q^ , the square of the 
4- momentum transfer between the K and the tt. This in- 
volves calculating a so-called 3-point function, shown in 
Figure [II where two different hadron operators are pulled 
apart in time and a current operator is inserted to con- 
vert a quark from one flavor to another. This is a much 
more complicated calculation than that of a 2-point func- 
tion and in addition must be done at a range of spatial 
momenta so that the q^ dependence of the form factor 
can be extracted. The result required for A' ^ tt is the 
form factor at g^ = and this can be obtained cither by 
extrapolation from results at the lattice spatial momenta 
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available or by generating appropriate spatial momenta 
to give g^ = by creative use of the lattice boundary 
conditions. The RBC/UKQCD collaboration have new 
results this year on the K -^ -k form factor using these 
techniques and the domain wall quark formalism 14 1. 
The advantage of using K semileptonic decay for Vus is 
that the chiral extrapolation of the form factor is known 
to be relatively benign because of the AdemoUo-Gatto 
theorem. This has been used in the past to estimate the 
difference of the form factor from 1.0 at q^ = (and 
this is what existing Vus determinations are based on) 
but it is to be hoped that lattice QCD can give a more 
accurate result. The RBC/UKQCD collaboration find 
/+(0) = 0.9644(49) which yields K, = 0.2249(14) ^. 
The error in the lattice result is estimated to be 0.6% 
with results available currently at one value of the lat- 
tice spacing. Results from multiple values of the lat- 
tice spacing should allow this to be improved. Fig- 
ure [8] shows their results for /+ (0) as a function of u/d 
quark mass (given by mj) and their chiral extrapolation 
compared to continuum estimates from different groups. 
Both the K12 and Kli results show that lattice QCD 
has an important role to play in the determination of 
Vus- This is summarised in the Flavianet plot [29|; Fig- 
ure [9l that shows the "K, ., from Kl'i (using lattice results 
from RBC/UKQCD [3]) and the ratio of Vus/Vud (us- 
ing HPQCD lattice results [13] ) compared to first-row 
unitarity constraints. 




FIG. 7: A 3-point function in lattice QCD involves two 
hadrons separated in time with the insertion of an operator in 
between to convert a quark from one flavor to another. Fit- 
ting this correlator as a function of the time position of the 
current, and of the time separation of the two hadrons, for var- 
ious values of the spatial momenta of the two hadrons then 
allows the form factor as a function of squared 4-momentum 
transfer, (f' , to be determined. 



Neutral kaon mixing is an important process through 
which we can study CP violation in that system. In the 
Standard Model it proceeds via the 'box' diagram shown 
in Figure [To] with W boson exchange. On the lattice, we 
are working at relatively low energy scales compared to 
the W boson mass and the version of the box diagram 
that is appropriate is that of a 4-quark operator from the 
effective weak Hamiltonian at low energy scales. The ma- 
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and the two numbers quoted separately. Jk on its own is 
a relatively simple calculation, as described above, and 
has been calculated accurately in lattice QCD by sev- 
eral groups. The 4-q operator is a harder calculation and 
is complicated to renormalise in formalisms that do not 
have the complete continuum chiral symmetry. This has 
made it a good calculation for the domain wall quark 
formalism. 





FIG. 8: The scalar form factor, /o(0), for K —> tt decay 
obtained from lattic e Q CD calculations with 2+1 flavors of 
domain wall quarks [2^. The plot shows the form factor as 
a function of quark mass and fits used to extrapolate to the 
physical point. For comparison values are given from contin- 
uum estimates of various groups. 
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FIG. 9: A combined fit by the Flavianet WG to test con- 
sistency of Vus and Vus/Vud determinations using K13 and 
K12 decays along with lattice results. The constraint from 
first-row unitarity of the CKM matrix is also shown. 



trix element of this operator between neutral kaon states 
is parameterised by f^BK where fx is the decay con- 
stant and Bk, the amount by which the matrix element 
differs from /|- is known as the bag factor. An evaluation 
of the 4-q operator matrix element in lattice QCD along 
with experimental results for the direct CP violation ra- 
tio of 27r decay rates of Kl and Ks, e, allows a constraint 
on the CKM elements that enter the 4 — g operator from 
c and t quarks appearing in the box diagram [30(] . There 
is no need in principle to separate the matrix element 
into its fn and Bk components, but this is usually done 



FIG. 10: The 'box diagram' reduces to a 4-quark operator 
in the effective weak Hamiltonian. It is the matrix element of 
this latter operator that is calculated on the lattice. 

RBC/UKQCD have given a result for Bk using do- 
main wall quarks this year on their configurations in- 
cluding M, d, and s sea quarks. They have calculations 
on both 16'^ and 24^^ lattices at a single lattice spacing 
of O.llfm. They obtain Bf^(2GeV) = 0.524(10)(28) 
where the first error is statistical and the second sys- 
tematic [6|. However, their result for Jk on the same 
lattices is currently 5% low. Encouraging preliminary re- 
sults, but no final numbers, were also presented this year 
from a mixed-action approach using domain wall valence 
quarks on the MILC gluon configurations that include sea 
quarks with the improved staggered formalism at multi- 
ple values of the lattice spacing [31[. Errors below 5% 
are clearly possible on Bk from lattice calculations in 
the near future. 

Charm physics provides a key test of lattice QCD cal- 
culations because, as described earlier, it is more sensi- 
tive to discretisation errors coming from the lattice QCD 
method than light quark physics. A new formalism called 
HISQ quarks [5| has shown this year that it is possible to 
treat charm quarks in essentially the same way as light 
quarks. This makes for a very fast method (because it 
is based on the fast staggered formalism, but working 
at larger quark masses it is even faster) with the added 
advantage of being able to make use of light quark chi- 
ral symmetry to give conserved currents for calculating 
decay rates that do not need renormalisation (and their 
associated systematic errors). This then gives errors at 
the few percent level [12| . 

One significant test of the formalism, and one that has 
not been available from previous lattice calculations, is 
that of the simultaneous determination of the spectrum 
of charmonium states and charm-light states with the 
same charm quark propagators. Of course, in QCD we 
know that there is only one charm quark with a given 
mass, but most approximations to QCD, such as po- 
tential models, find it impossible to handle both sets of 



states in the same approximation since their internal dy- 
namics is very different. On the lattice it is also true 
that systematic errors are very different in the two sys- 
tems, with charmonium states being more sensitive to 
discretisation errors, so that a very accurate discretisa- 
tion is needed to be able to describe both successfully. 
With HISQ we have this and so the value that the charm 
quark mass needs to take is fixed to get the rjc mass cor- 
rect (this being the lowest-lying charmonium state and 
one whose mass is most accurately calculated on the lat- 
tice). The Dd and Ds masses are then non-trivial pre- 
dictions given a u/d mass from ttItt and an s quark mass 
from TTiK- Very accurate results are obtained using va- 
lence HISQ quarks on the MILC ensembles, which agree 
with experiment with 6 MeV errors, see Figure [TlJ This 
level of accuracy requires an understanding of corrections 
to the meson masses in the real world from QED effects 
and the fact that the u and d masses are not the same. 
To be working at the level of precision where these effects 
have to be considered is very exciting. 
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FIG. 11; Masses of the D'^ and Ds meson as a function of 
the u/d mass in units of the s mass at three values of the 
lattice spacing using the MILC ensembles. The lines gives 
the simultaneous chiral fits to all three lattic e sp acings and 
the dashed line the continuum extrapolation [I^l. The final 
error bars are given by the shaded bands, offset from the 
dashed lines by an estimate of electromagnetic, rriu 7^ nid and 
other systematic corrections to the masses. The experimental 
results are marked at the physical nid/rns- 



The D and Dg decay constants can also be determined 
to an accuracy of 2% using exactly the same method as 
for fx and /^ described above. The results are shown in 
Figure fT2l The leptonic decay rates of D and Dg mesons 
have now been measured by experiment and, given a 
value for the appropriate CKM element from elsewhere, 
can be converted into a value for the decay constant that 
can be compared to lattice results. Figure [12] compares 



such experimental determinations of fn^ from BaBar, 
Belle and CLEO-c, using Vcs = Vud, with lattice results. 
The Fermilab/MILC collaboration lattice result is also 
shown. They use the clover formalism for charm quarks 
on the MILC ensembles and produced first results for fo 
and /ds ahead of experimental results and with errors of 
7% [32| . The value of fo^ shown in figure [13] is an up- 
dated one from this year's lattice conference of 254(14) 
MeV [33|. The calculation of fu and fo^ in the clover 
formalism does require a renormalisation (and the error 
associated with this is included in the error estimate) 
and the mass of the charm quark in this case is fixed 
from the Dg mass itself, so further independent checks of 
this calculation against other quantities need to be done. 
The two lattice calculations agree well although they use 
very different charm quark formalisms, albeit on the same 
gluon ensembles. Further work is underway from other 
groups using other formalisms and ensembles. 

On this quantity the lattice results arc ahead of ex- 
periment, although experimental errors arc expected to 
improve by a factor of two over the coming year. The 
experimental central values will have to come down as 
that happens if there is to be agreement. One outstand- 
ing issue is that of electromagnetic effects in the experi- 
mental result since the decay constant is defined in pure 
QCD [3J|. More work needs to be done in the D and Dg 
cases to ensure that is not an issue at the level of preci- 
sion we are aiming for in this comparison between theory 
and experiment. 

Neutral B mesons can mix in an analogous way to 
K mesons. In this case it is simpler, being dominated 
by a single box diagram with t quarks in, and reduc- 
ing to a 4-quark operator with coefficient Vtdjs^th (^^^ 
Figure fTO|) . The matrix clement of this operator is pa- 
rameterised by /^Bb where /b is the decay constant. 
In the past there has been a lot more work on the de- 
cay constant than on the 4-quark operator. This year, 
however, two groups, the HPQCD collaboration and the 
Fermilab/MILC collaboration, presented new results for 
the ratio of the 4-quark operator matrix element for the 
Bg to that of the Bd ;35l| . The results are at a preliminary 
stage currently but will be improved significantly over the 
The ratio ^ 



commg year, ine ratio ^ = /b, \/BbJ fB\/BB is more 
accurate than the individual numbers because we cannot 
escape a renormalisation (and its associated error) of the 
matrix elements in this case, but the renormalisation is 
the same for B and Bg. Previously the HPQCD collabo- 
ration gave a result for /bs/Jb oi 1.20(3) i36t]. This year 
the Fermilab/MILC collaboration give their value for this 
ratio as 1.22(3) [33[. Our prejudice is that ^ should be 
very similar to /b^ / fs, and the aim is certainly to achieve 
the same level of error, or better. HPQCD have given a 



resuh for fB,\J Bb, of 0.281(21) GcV [371, where a lot of 
this error comes from the renormalisation. 

Exclusive B semileptonic decay is an important route 
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FIG. 12: Results for the D, Ds, K and tt decay constants 
on the very coarse, coarse and fine MILC ensembles using 
valence HISQ quarks. The chiral fits are performed simulta- 
neously with those of the corresponding meson masses and 
the resulting continuum extrapolation curve is given by the 
dashed line. The shaded band is the final result. Experimen- 
tal results are shown on the left. For K and vr these are from 
the PDG [131. For L» and Ds these are from CLEO-c J^, Hi] 
(with ^ and r results shown separately) and BaBar [ioj ]. 
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FIG. 13: A comparison of lattice results for fos from calcula- 
tions including 2-f 1 fiavors of sea quarks, compared to experi- 
mental determinations from the leptonic decay rate and using 
Vcs ~ Vud- The two CLEO-c results from /i and r channels 
are separated. 
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to determining the CKM elements Vcb and V^b- This 
year the Fermilab/MILC collaboration have determined 
the appropriate form factor for the process B — > D*^, in- 
cluding the effect of sea quarks fuUy for the first time [43] . 
The rate of this decay at zero recoil is proprotional to the 
square of Vcb x ^Ai(l), where /iai(1) is the form factor 
for B and D* relatively at rest. Figure [U shows the Fer- 
milab/MILC results for ft,^j(l) at three different values 
of the lattice spacing using the MILC gluon field config- 
urations including the effect of sea quarks. From their 
results they determine /iai(1) = 0.942(12)(19) where the 
first error is statistical and the second systematic. Us- 
ing the HFAG experimental average [43[, this leads to a 
value for Vcb of 38.7(0.7) (0.9) x lO'^ where the first error 
is from experiment and the second from the lattice. 

Vub can be extracted from the exclusive B -^ ttIv pro- 
cess. Here lattice results are still not very accurate be- 
cause the important kinematic region is one in which 
the n mesons arc moving rather fast and this gives a 
very noisy signal on the lattice. Work is underway to 
ameliorate this 44|. A recent theoretical determination 



FIG. 14: The form factor at zero recoil for B -^ D*lv cal- 
culated by the Fermilab/MILC collaboration on the MILC 
configurations at three values of the lattice spacing (in fm^ . 
Good consistency between them is observed. The final result 
is given as 0.924(12) (19) HI]. 



of Vub using combined lattice results from HPQCD and 
Fermilab/MILC along with light-cone sum rules gives 
Vub = 3.47(29) X 10-3 ^. The issue of compatibility of 
this result with that from the inclusive b —^ u decay is 
becoming an important one [46[. 



CONCLUSIONS 

Lattice calculations including the full effect oi u, d 
and s sea quarks are in excellent shape. Calculations 
using improved staggered quarks continue to get better 
and new results are now appearing from other quark for- 
malisms. There have been significant new results this 
year in strange and charm physics. I have concentrated 
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on results relevant to CKM physics and have not men- 
tioned many other areas of progress in lattice calcula- 
tions. The reader is referred for those to the proceedings 
of this year's lattice conference [3l- 

In Figure [T5l I have collected recent lattice results into 
a set of constraints on the upper vertex of the standard 
unitarity triangle plot. The lattice inputs needed are Bk, 
Jk/U, f+{K ^ TTli^), F{B ^ D*lv), f+{B ^ TTlu) and 
fB^\/BBj Ib^/Bb- It is important to use lattice results 
from the calculations including the full effect of u, d and 
s sea quarks. Old results in the quenched approximation 
are not reliable enough for this. They have in any case 
now been superseded and should not be used. 

In the next two years the lattice errors on these CKM 
constraints should halve. The robustness of our error 
estimates will be further tested against experiment using 
other gold-plated results. Fg+g- for ip and T are good 
tests for c and h physics, for example. We have no free 
parameters when we do this and so it is a very stringent 
test. The era of precision lattice QCD calculations has 
really arrived. 



de Water and Kit Wong. 
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FIG. 15: Constraints on the unitarity triangle from the CKM 
matrix using current lattice results that include the effect of 
u, d and s sea quarks. The lattice errors dominate these 
constraints and they can be halved over the next two years. 
The cross gives the current constraint on the vertex quoted 
in the PDG [13]. 
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